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Intensity  vectors  of  the  Pekeris  a  -0*1  are  calculated  from  values  of 
acoustic  pressure  and  particle  velocity  which  are  obtained  by  numerical 
integration.  Plots  of  the  vectors,  for  the  cases  of  no-seabed,  a  rigid  seabed 
and  fast  and  slow  seabeds.  Illustrate  the  potential  usefulness  of  the 
intensity  vector  approach  as  an  aid  to  understanding  the  physics  of  underwater 
sound  propagation. 
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INTRODUCTION 


Sound  propagation  in  the  saa  can  ba  a  coaplax  procaaa  Cor  which 
mathematical  aodals  of  varying  complexity  hava  bean  developed.  The  Pekarla 
nodal  ia  the  noat  basic  o£  these;  it  consists  of  a  hoaogeneous  layer  of  fluid, 
which  contains  a  tine  harmonic  point  source  of  sound,  lying  over  an  infinite 
half-space  of  another  hoaogeneous  fluid  [1-3].  Numerical  values  of  pressure 
obtained  from  this  nodal  are  usually  presented  as  plots  of  propagation- loss 
versus  range,  at  a  fixed  depth;  or  alternatively,  as  a  contour  plot  of 
transnission  loss  over  depth  and  range.  Jensen  £  Kupeman  [4]  use  both  these 
presentations  in  their  report  which  compares  nunerical  results  which  were 
obtained  fron  various  theoretical  propagation  node Is,  which  are  also 
described.  Little  progress  has  been  reported  in  the  presentation  of 
experinental  neasurensnts  since  Mood's  scale  nodal  work  [S]. 


Plots  of  acoustic  Intensity  vectors  Show  not  only  the  magnitude  of  the 
energy  flaw  but  also  its  direction.  Bence,  intensity  vector  plots  enable  an 
analyst  to  distinguish  between  cases  in  which  there  is  a  net  outflow  of 
energy,  and  those  cases  in  which  the  energy  circulates  in  closed  loops. 
Pettersen  [6]  has  presented  plots  which  show  energy  flowing  fron  sane  simple 
source  and  sink  configurations,  and  also  plots  which  show  diffraction  around 
and  reflection  fron  a  thin  barrier.  Spicer  [7]  and  Janes  [8]  have  computed 
acoustic  Intensity  vectors,  for  a  fluid-loaded  plate  and  shell,  respectively, 
which  help  to  illustrate  the  physics  of  fluid-structure  interaction. 
Publications  by  Fahy  [9]  and  Reinhart  a  Crocker  [10]  have  demonstrated  that 
the  'two-microphone  technique'  Is  now  a  practical  tool  for  the  experinental 
measurement  of  acoustic  intensity. 


In  this  memorandum,  the  mathematical  formulae  needed  to  evaluate 
numerically  the  pressure  and  velocity  fields,  of  the  Pekerls  model,  are  given 
for  the  special  cases  of  no-seabed,  a  rigid-seabed,  and  fast  and  slow  seabeds. 
The  calculation  of  intensities  from  numerical  values  of  acoustic  pressure  and 
particle  velocity  is  straightforward.  Plots  of  intensity  vectors  at  low 
frequencies  and  over  short-ranges  are  presented  and  discussed.  These 
intensity  vector  plots  reinforce  previous  numerical  work  on  the  Pekerls 
model,  and  it  is  hoped  that  they  will  prove  to  be  of  sufficient  Interest  to 
encourage  similar  presentations  in  work  on  more  realistic  models,  the 
Mathematics  necessary  to  extend  the  work  contained  here  to  the  model 
consisting  of  an  arbitrary  number  of  homogeneous  layers  of  fluids  and  elastic 
solids  is  described  elsewhere  [11,12],  In  the  more  realistic  models  a  rich 
spectrum  of  propagating  waves,  such  as  the  Stone ley  and  Rayleigh  Interface 
waves,  may  be  excited  and  it  is  expected  that  intensity  vector  plots  may 
reveal  which  of  these  propagating  waves  are  likely  to  be  detected 
experimentally. 


2.  PROBLEM  FORMULATION 


The  Hlwrli  Modal  for  sound  propagation  In  a  shallow  ocaan  consists 
of  a  tlas  harmonic  point  sourca  of  sound,  of  fraa-flald  prasaurs 
Po*xp(  i*1R0-iwt  )/Rq,  which  is  located  in  a  layer  of  homogeneous  acoustic 
fluid,  whose  upper  boundary  t-N  is  in  contact  with  a  mcuia,  and  whose 
lower  boundary  x-O  is  in  contact  with  a  half-space  of  a  second  homogeneous 
acoustic  fluid.  This  geometry  is  shown  in  Figure  1. 

The  differential  equations  satisfied  by  the  acoustic  pressure, 
expressed  in  cylindrical  polar  coordinates  and  omitting  the  time  variation 
exp(-iwt),  for  this  axisyomtrlc  problem  axe 


v2p1<r,s)  +  1t2p1<r,s)  -  -4w8(r)a(*-*0)/r  (2.1) 

V2P2<*.*)  +  *^P2(r.*)  -  0  (2.2) 

The  boundary  conditions  of  continuity  of  pressure  at  the  boundaries  of 
the  fluid  layer  are 


p^r.B)  -  0  (2.3) 

Pj( r,0)  -  P2(r,0)  (2.4) 

and  continuity  of  normal  displacement  at  the  interface  between  the  fluid 
layer  and  the  fluid  half-space  requires  that 


( VP^CPp^r.zya*]  -  (l/P2)(»P2(r,  *)/«*]  (2.5) 

s*o  S-O 


The  remaining  boundary  conditions  are  those  Ispossd  by  the  radiation 
condition,  vis.,  at  large  distances  from  the  source  the  pressure  consists 
only  of  outgoing  waves.  The  solutions  of  equations  (2.1)  -  (2.5)  axe 
well  known  and  are  swmarised  in  the  Appendix  for  special  cases  of  the  seabed 
constants . 


The  time  averaged  acoustic  intensity  vectors  in  the  fluid  layer  have 
radial  and  vertical  components  defined  by 


Ilr(r,s)  -  ( 1/2 )  Re[pl(r,s)Vlr(r,s))  (2.6) 

Ils(r,*)  "  (V2) 


(2.7) 


The  amplitude  and  phase  of  the  vectors  are 


*  tr?r(r'*>  +  I!s(e'*)]1/2  (2l9) 

&t(r,z)  -  tan-1tiia(r,s)/Itr(r,*)]  (2.9) 

Equations  (2.6)  -  (2.9)  also  apply  to  the  fluid  half-space  If  the 
subscript  1  Is  replaced  by  the  subscript  2. 


3 .  NUMERICAL  RESULTS 


(a)  General 


Fortran  computer  prograss  have  been  written  to  calculate  and  store 
on  a  disk  file  the  numerical  values  of  pressure  and  particle  velocities 
at  points  on  a  selected  grid  In  the  r-s  plane,  for  the  various  types  of 
seabed  considered  here.  A  separate  Fortran  program  uses  this  data  to  fora 
and  plot  the  time  averaged  intensity  vectors. 


The  numerical  constants  In  SI  units  used  to  produce  the  results  shown 
In  Figures  2-13  area 


Fluid  Layer 


pt  -1000.0 


C j  -1500.0 


*0  “55.0 
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Fast  Fluid  Seabed  p2  -1500.0 


c2  -1650.0 


Slow  Fluid  Seabed  p2  -1500.0 


c2  -1350.0 


The  Integrals  In  equations  (A7-A15)  must  be  evaluated  numerically 
because  cloeed  form  expressions  are  not  available.  Because  the  Hankel 
transform  representation  Is  arbitrary  with  respect  to  a  solution  of  the 
homogeneous  wave  equation,  it  is  necessary  to  introduce  damping  Into  the 
system  In  order  to  lift  the  poles  off  the  real-axis.  Damping  Is  Included 
by  setting  the  sound  speeds  In  the  fluids  to  the  complex  values 


cl  *  cl<l-lTll> 

e2  *  c2(1-1tV 


(3.1) 


where  t)j 


and  r\2  are  Chosen  as  0.001.  This  value  is  sufficiently  small 


/ 


to  avoid  significant  attenuation  in  the  fluids  over  the  distances  considered, 
yet  it  is  large  enough  to  slice  numerical  integration  without  ill-conditioning 
being  present,  the  masexical  approximation  of  the  Integrals,  after  truncating 
the  infinite  limit  to  a  sufficiently  high  value,  is  accomplished  by  a  simple 
adaptive  Gaussian  quadrature  schema .  Parallel  computation  on  an  array  of 
r-valuee  results  in  a  considerable  saving  in  coeputar  time. 


In  Figures  4-13  the  plotted  lengths  of  the  normalised  vectors  are 
proportional  to  RgV(l),  Which  is  a  spherical  spreading  correction. 
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Figure  2(a)  shews  the  wavenumber-frequency  plot  of  the  rigid  aaabad 
where  fraa  wavw  cut-on  at  frequencies  4.7Ha,  14. IBs  and  23 . 4Hx .  Figure  2(b) 
•hows  tha  wavenustoer- frequency  plot  of  tha  fast  aaabad,  where  a  aingla  branch 
cut a— on  at  11 . 3  Hz .  Tha  linaa  and  o“k2 »  batvaan  which  tha  raal 

wavanuabars  for  tha  fast  aaabad  ara  conatrainad  to  lia,  ara  included  on  tha 
plot  as  reference  linaa. 


(c)  Tha  Reflection  Coefficient 


Tha  plane -wave  reflection  coefficient  at  the  seabed 


R<®)  -  CpjCjCosS-p^Cj-c^ain^n/tp^eoae+PjVtCj-cjjain2*)]  (3.4) 

is  of  interest  because  the  transform  representation  of  the  pressure  field 
resolves  the  pressure  in  tha  fluid  layer  into  a  continuous  spectrum  of 
waves  propagating  both  upwards  and  downwards  at  angles  of  incidence  e  given 
by  the  equation 


a  -  ktsine 


In  Figure  3(a)  the  magnitude  of  the  reflection  coefficient  for  tha 
slow  seabed  decreases  steadily  with  increasing  8  until  it  vanishes  at 
the  angle  of  perfect  transmission,  8-64°.  At  larger  values  of  8  the 
magnitude  of  the  reflection  coefficient  increases  steadily,  attaining  its 
ma-rljssi  value  of  unity  at  8—90° . 

In  Figure  3(b)  the  magnitude  of  the  reflection  coefficient  for  the 
fast  seabed  is  an  almost  constant  o.2S  until  0-40° .  it  then  increases 
rapidly  to  the  value  of  unity  at  8-65°,  tha  case  of  perfect  reflection, 
after  which  it  remains  at  unity. 


(d)  Intensity  Vector  Plots  at  2.5Hz 


Figures  4  and  5  shew  tha  intensity  vector  plots  at  a  frequency  of 
2.5Hz  where  tha  wavenumber- frequency  plots  of  Figure  2  show  that  there 
are  no  propagating  modes.  There  1s  little  difference  between  the  plots 
for  no  seabed,  tha  slow  seabed  and  tha  fast  seabed,  all  of  Which  show  the 
characteristics  of  a  dipole  sound  source  [6]  where  most  of  the  sound  energy 
is  transmitted  into  tha  seabed.  Tha  case  of  the  rigid  seabed  is  somewhat 
special  because  tha  intensity  is  only  non-zero  by  virtue  of  tha  finite  value 
of  m •  Hare  tha  energy,  which  is  two  orders  of  magnitude  smaller  than  that 
of  tha  other  plots,  is  being  absorbed  in  tha  layer. 
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(e)  Intensity  Vector  Plots  at  s.OBs 


Figures  6  and  7  show  the  Intensity  vector  plots  at  a  frequency  of 
5.0Hz,  where  the  wavenumber-frequency  plots  of  Figure  2  show  that  a 
propagating  node  is  present  only  in  the  rigid  seabed  case.  This  node 
quickly  becomes  dominant  as  r  Increases,  its  Intensity  vector  components 
being 


I  ~  cos2(»rz/2H)/r,  I  -  o 

E  Z 


Again,  there  is  little  difference  between  the  plots  for  no  seabed,  the  slew 
seabed  and  the  fast  seabed,  where  the  Intensity  vectors  exhibit  the 
Characteristics  of  dipole  radiation  of  energy  into  the  seabed. 


(  f )  intensity  Vector  Plots  at  15 . OKz 


Figures  8  and  9  show  the  intensity  vector  plots  at  a  frequency  of 
15.0Hz,  where  the  wavenusber-frequency  plots  of  Figure  2  show  that  two 
propagating  nodes  are  present  In  the  rigid  seabed  case,  and  one  is 
present  in  the  fast  seabed  case.  Again  the  plots  for  no  seabed  and  the 
slow  seabed  approximate  to  dipole  radiation  of  energy  into  the  seabed. 


In  Figure  8(b)  the  two  propagating  modes  coeblne  to  produce  an 
interference  pattern  which  shows  energy  propagating  along  the  layer  with 
non-zero  vertical  intensity  vector  components  although  the  intensity  vectors 
of  the  separate  modes  are  everywhere  horizontal.  This  Is  because  intensities 
do  not  add  linearly  -  the  pressure  of  the  first  mode  Interacts  with  the 
velocity  of  the  second  mode  and  vice  versa,  Schultz  (131.  The  total  energy 
flow  in,  however,  the  sum  of  the  energies  contained  In  the  individual  modes 
because  the  contribution  of  the  cross-terms,  when  integrated  across  the  layer 
thickness.  Is  zero. 


In  Figure  9(b),  at  small  distances  from  the  source  the  intensity  vectors 
approximate  to  those  of  dipole  radiation  of  energy  Into  the  seabed.  At 
larger  distances  the  propagating  mode,  whose  cut-on  frequency  is  11.3Hz, 
interacts  with  an  evanescent  mode  to  produce  a  closed  loop  of  energy 
circulation.  The  propagating  mode  may  be  resolved  Into  two  waves,  one  of 
which  propagates  upward  and  the  other  downward,  their  angle  of  incidence  being 
greater  than  the  65*  at  which  perflect  reflection  occurs.  This  mode  just 
begins  to  dominate  at  the  maximum  range  shown,  viz.  150m. 
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( g )  intensity  Vector  Plots  at  25.0Hz 


Figures  10  and  11  show  the  intensity  vector  plots  at  a  frequency  of 
25.0Hz  where  the  wavenumber- frequency  plots  of  Figure  2  show  that  three 
propagating  nodes  are  present  in  the  rigid  seabed  case  and  one  is  present 
in  the  fast  seabed  case.  The  plots  for  no  seabed  and  the  slow  seabed  are 
similar  and  are  typical  of  intensity  patterns  obtained  froa  well-separated 
sources. 


In  Figure  10(b)  the  three  propagating  nodes  conbine  to  produce  an 
interference  pattern  which  shows  energy  propagating  horizontally  at  the 
top  and  bottan  of  the  layer.  The  plot  for  the  fast  seabed.  Figure  11(b), 
is  not  unlike  the  plots  of  Figure  10(a)  and  11(a);  the  single  propagating 
node  does  not  doninate  the  solution  until  radial  distances  are  In  excess 
of  300n  (not  shown). 


(h)  Intensity  Vector  Plots  in  the  Seabed 


Figures  12  and  13  show  the  intensity  vector  plots,  at  a  frequency  of 
25.0Hz,  in  the  slow  and  fast  seabads  respectively.  The  reduction  in  Magnitude 
of  the  vectors  as  the  interface  is  crossed  is  presumably  due  to  partial 
reflection  at  the  seabed.  As  predicted  by  Snell's  law,  the  intensity  vectors 
in  the  slow  seabed  are  refracted  slightly  towards  the  normal,  while  the 
Intensity  vectors  in  the  fast  seabed  are  refracted  away  from  the  normal. 


4.  CONCLUDING  REMARKS 


Formulae  have  been  given  from  which  values  of  the  acoustic  intensity 
vectors  of  the  Pekerls  model  of  underwater  sound  propagation  may  be  calculated 
numerically.  The  limited  number  of  plots  shown  here  illustrate  the  way  in 
which  both  the  magnitude  and  direction  of  the  energy  flow  in  the  sea  can  be 
controlled  by  the  properties  of  the  seabed,  especially  when  these  properties 
are  such  as  to  allow  the  propagation  of  normal  modes.  Of  particular  interest 
are  the  complex  interference  patterns  and  circulatory  energy  flows  which  this 
new  method  of  presenting  propagation  data  highlights.  These  preliminary 
intensity  vector  plots  are  considered  to  be  of  sufficient  interest  to  justify 
further  numerical  work  using  more  realistic  models  of  the  seabed  in  which, 
for  example,  shear  wave  motion  is  included.  The  practical  value  of  intensity, 
as  opposed  to  pressure,  measurement  could  be  investigated  by  means  of 
interacting  theoretical  and  scale  model  work. 


E  A  Skelton  (  hso) 


REFERENCES 


1.  BREKBOVSKIKH  L.M. ,  waves  in  Layered  Media,  Academic  Press,  1980. 

2.  BREKBOVSKIKH  L. ,  lysanov  Yu. ,  Fundamentals  Of  ocean  Acoustics, 

Spr inger-Ver lag ,  1982. 

3.  GABRIELSON  T.B.,  Mathematical  Foundations  for  Normal  Mode  Modelling 
in  Waveguides ,  Naval  Air  Development  centre,  Pennsylvania, 

Report  No.  NADC-8 1284-30,  1980. 

4.  jensen  F.B.,  KUPERMAN  w.A. ,  Consistency  Tests  of  Acoustic  Propagation 
Models,  saclant  Research  Centre,  La  spezia,  Italy,  Memorandum  SM-157, 
March  1982. 

5.  NOOO  A.B.,  Model  Experiments  on  Sound  Propagation  in  Shallow  Seas, 

J.  Acoust.  Soc.  Am.,  31(9),  ppl2l3— 1235 ,  September  1959. 

6.  PETTERSEN  O. ,  Sound  Intensity  Measurements  for  describing  Acoustic 
Power  Flow,  Applied  Acoustics,  14,  pp387-397,  1981. 

7.  SPICER  W.J.,  Acoustic  Intensity  Vectors  from  an  Infinite  Plate  with 
Line  Attachments,  Admiralty  Marine  Technology  Establishment,  Teddington, 
AMTE(N)TM81086,  October  1981. 

8.  JAMES  J.H.,  Acoustic  Intensity  Vectors  inside  Infinite  Cylindrical 
Elastic  Shell,  Admiralty  Karine  Technology  Establishment ,  Teddington, 
AMTE(N)TM82103,  December  1982. 

9.  FAHY  F.J.,  Measurements  with  an  Intensity  Meter  of  the  Acoustic  Power 
of  a  Small  Machine  in  a  Room,  J.  Sound.  Vib.,  57(3),  pp311-322,  1978. 

10.  REINHART  T.E.,  CROCKER  M.J.,  Source  Identification  on  a  Diesel  Engine 
using  Acoustic  Intensity  Measurements,  Noise  Control  Engineering, 

18(3),  pp84— 92,  1982. 

11.  SPICER  W.J.,  Free-wave  Propagation  in  and  Sound  Radiation  by  Layered 
Media  with  Flow,  Admiralty  Karine  Technology  Establishment,  Teddington, 
AKTE(N)TM82102,  December  1982. 

12.  SKELTON  E.A.,  Sound  Radiation  from  a  Cylindrical  Pipe  Composed  of 
Concentric  Layers  of  Fluids  and  Elastic  Solids,  Admiralty  Marine 
Technology  Establishment,  Teddington,  AtrrE(N)TM83007,  January  1983. 

13.  SCHULTZ  T.J.,  et  al.  Measurement  of  Acoustic  Intensity  in  a  Reactive 
Sound  Field,  J.  Acoust.  Soc.  Am.,  57(6)  Part  1,  pp!263-1268,  1975. 


-15- 


FREQUENCY (H* ) 


RADIAL  WAVENUMBER  v.  FREQUENCY  PLOTS 


•U- 


ANTE  (N  )  TM*J0  M 


(a)  No  Saabad 


-•  -«  ..  «  ,  ,  ,  ,  ,  ,  ,  ,  . 

■*  •«  •,  i  -«  •>  ■■  ..  ,  ,  ,  ,  •  ,  ,  , 

-•  •»  -•  ■»  ,  <  •>  ■.  ,  , 

X  ir  \  1  «»  -• 


(b)  Rig  Id  Sea bad 


/»-***■«■* 


'■’•'»  ■•  -  -  -  -  *  * 

'»'»'*■'  •*  >•  -•  -  -  -  - 

-.  -  -  -  - 

•*  ■•  ■•  •• 
-  -  - 
"•  *•  -  - 

■i  \  -  - 

—  \  -•  ■•  -•  ■• 


INTENSITY  VECTORS  AT  2  5HI 


f 

4 


-a o- 


amtb <n) TNaaoes 


(a)  Slow  Seabed 


■•••■>  t  <  •  ,  -  .  ••  .  .  -  -  .  .  . 

-•'»  i  ■»  -•  ■.  ■■  ,  •  - 

-»  ■»  ■»  ■»  ■»  •»  -•  ■»  -»  -  ->  - 
-»  ->  •»  -  -  - 

••  ■»  •»  ■*  •» 

V  V  W'i'n'llSNN'*'*  'H't'l  S'!  1  ->»  -» 

X  >*  >»  -»  -» 

X  X  ^  -» 

X  X  i  V  « 

'=0  r*188 


(b)  Fast  Seabed 


'SV*'S'SN»  N  N  "  ~  '  -  *  '  *  *  ' 

"*  ■»  •*  ■»  ”•  ■•  "*  "»  ■*  '•  '  "*  ' 

’  '»'*  -  ■*  -  -  **  '  ■* 

"*  -*  "*  "*  "•  "•  "• 

XXXXXX^'b'b'bN*v»'b''»'b>»'»'»>»'»^'»'*^>»>»^>»>» 

r*188 


FIG.  5  INTENSITY  VECTORS  AT  2.  5Ha 
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FIG.  7  INTENSITY  VECTORS  AT  5Hz 
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FIG.  8  INTENSITY  VECTORS  AT  ?5Hz 
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FIQ.  9  INTENSITY  VECTORS  AT  15Hz 
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FIG.  11  INTENSITY  VECTORS  AT  25H* 
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FIG.  12  INTENSITY  VECTORS  AT  25Hz 
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FIQ.  13  INTENSITY  VECTORS  AT  25Hz 


amti <n> TMeeoes 


APPENDIX 


Formula*  for  acoustic  pressure  and  partlc la  velocity  of  the  Pekerls  nodal 


(a)  No  Seabed 


In  the  special  case  in  which  the  layer  and  the  half-space  consist  of 
Identical  fluids,  the  pressure  In  the  fluid  Is  obtained  by  the  method  of 


Pjtr.x)  -  pQe*p(lX1R0)/R0  -  pQe*p(  iX^  )/Rt 


(AX) 


where 


2  2  2 
R0  *  *  +  <*'V 

2  2  2 
R*  -  r  +  <*+*0-2H) 


The  radial  and  vertical  particle  velocities  are 

vir(r'*)  "  (P0*/i“P1X(XJe)l-l/R0)«>5><Ut,l«())/»0  -  (lkj-l/Rj  )**P<iX1R1  )/R*}  (A2) 


Vljt(  r,a)  -  (Pg/x «p4  ) {( *-*0  )( lXj  -1/R0  )*ap(  1XJR0  )/R* 


-  (s+Sjj-XHXUCj-l/RjJeapdk^)/**)  (A3) 


(b)  Rigid  Seabed 


The  rigid  seabed  nay  be  considered  as  the  limiting  case  of  the  fluid 
half-space  as  p2  -  »  and  c2  -  •  .  the  well  known  series  solution  for  the 
pressure  In  the  fluid  layer  la  [1] 


p^r.aWfinp^H)^  cos(  t2n+l]Tfs0/2H)coa<  [2n+i)tr*/2H)H0(anr)  <A«) 

n-0 
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from  which  tha  radial  and  vertical  particle  velocitiee  nay  be  deduced  aa 


Vir<r,z)-(-2,rp0/<oPlH)  £  aftcoe(  C2n+l]rr*0/2H)coe(  [2n+l]rrz/2H)H1(anr )  (A5) 


Vlz(r,*H  -w  Pg/uiPjH  )  £  ( 2n+l  )coe(  [  2n+l  Irrx^ll  )ein(  [  2n+l  ]rrz/2H  )HQ(  «nr )  ( AS ) 


(c)  Past  or  Slow  Seabed 


In  the  general  case  of  tha  faet  seabed  (c2>Cj )  or  the  slow  seabed 
( c2<Cj ),  the  pressure  and  particle  velocities  in  the  fluid  layer  are  given 
in  the  region  0<z<z0  by 


P1<r,z)-2po  adt  D  Jo(car)sin(d1[B-*o]){/Jicos(diz)-ib02sin(01z)}do  (A7) 


v^Cr.iHZipg/up 


*  f  2-1-1 

1>  J  “  Pl  °  J1 


{  ar  )sin(  dt  [a-zQ]  )  {/Jlcos<  Pf  )-tbd2sin(  dj* )  }da 


Vla(r.z)-(2ip0/up1)  aO-1J0(enr)sin(S1[B-z(J]  HSlsin(0lz)+lb02cos(S1z)}da 


where  the  dispersion  relation  0  is 


D(a,w)-01cos(  d^-ibd^l^  djH) 


and  in  the  region  zQ<z<H  by 


|  o^V1. 


P,<r,tf2p0  afll  0  Jo(ar)sin(0llH-z]){Slcos<S1zo)-ib02sln(01zo))da  (A10) 


vlr(  r,  z )-( iXPg/up 


i>  j 


( ar  )sin<  0 1  [  H-Z  ] )  { d t cos(  <3 1  tQ  )-lbd2sin(  0 1  zQ )  }da 


-32' 


JQ<  or  >coa(  0^8-z]  ){01co«(0l*o  )-lb02«tn(0l* 


Vie(r.*)-(2ip0/«p1)  |  a nT1. 

0 

The  petuun  and  particle  velocities  in  the  lower  half-space  are 


P2(r 


.*)  ”  2P0  |  aO”1 


JQ(  otr  )ein(  0 1  [  8-*0  ]  )exp< -id2*  )do 


V2r(r.*)-(2ip0/up2)  |  a2D-lJ1(ar)ein(0i[H-*o]  )e*p(-i02z)da 


J  *»2d_1 


VM(r,z)  -  (— 2pQ/up2 )  I  «02D  *Jo(ar)sin(01[H-zo])e*p(-i02z)da 


)}da 

CM2) 

(MS) 

<M4) 

(MS) 
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